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Abstract 

We study gauge and gravitational field theories in which the gauge fix- 
ing conditions are imposed as constraints on classical fields. Quantization of 
fluctuations can be performed in a BRST invariant manner, while the main 
novelty is that the classical equations of motion admit solutions that are not 
present in the standard approach. Although the new solutions exist for both 
gauge and gravitational fields, one interesting example we consider in detail 
is constrained gravity endowed with a nonzero cosmological constant. This 
theory, unlike General Relativity, admits two maximally symmetric solutions 
one of which is a flat space, and another one is a curved-space solution of 
GR. We argue that, due to BRST symmetry, the classical solutions obtained 
in these theories are not ruined by quantum effects. We also comment on 
massive deformations of the constrained models. For both gauge and grav- 
ity fields we point out that the propagators of the massive quanta have soft 
ultraviolet behavior and smooth transition to the massless limit. However, 
nonlinear stability may require further modifications of the massive theories. 



1 Introduction and summary 



To quantize an electromagnetic field one could fix a gauge by imposing an operator 
constraint on physical states of the theory. For instance, in the Gupta-Bleuler 
approach one postulates 



where G(A) denotes a function of the gauge field operator A or its derivative, and 
is an arbitrary physical state of the theory. 

Alternatively, one could choose to impose the constraint on a gauge field already 
in a classical theory in such a way that ((H) is automatic upon quantization [H |2] . 
This can be achieved, e.g., by introducing into the classical action of the theory a 
Lagrange multiplier A times the function G(A) 



The field A has no kinetic or potential terms. Variation of the action w.r.t. A 
gives a constraint G(A(x)) = 0, which is just a classical counterpart of Small 
fluctuations of the fields in this theory can consistently be quantized for Abelian 
PI El 0] as well as non- Abelian and gravitational fields [I]. The resulting theories 
can be completed to be invariant under the Becchi-Rouet-Stora-Tyutin (BRTS) 
transformations 5J, and a Hilbert space of physical states can be defined by requiring 
that the states carry zero BRST and ghost charges Jlj. Quantum effects in the 
resulting theory are identical to those of the Gupta-Bleuler approach. 

Nevertheless, there is one difference in using the classically constrained theory 
that has not been explored. This difference could be seen in classical equations of 
motion. Variation of the action w.r.t. the gauge fields gives an equation of motion in 
which there are new terms proportional to A and/or its derivatives. Hence, classical 
field equations are modified, and, depending on allowed boundary conditions, new 
solutions could emerge. 

This becomes especially important for gravity. General Relativity (GR) with 
a non-zero cosmological constant does not admit Minkowski space as a solution 
of equations of motion. We will discuss classically constrained General Relativity 
(CGR) in section 3 and show that the latter does admit flat space as a solution even 
if the cosmological constant is not zero. This difference is clearly very important. 

Unimodular gravity (UGR) jH] is an interesting example of a partially constrained 
theory. In UGR the full reparametrization invariance of GR is restricted to a sub- 
group of volume preserving transformations. One practical difference between UGR 
and GR is that the cosmological constant problem in UGR is somewhat relaxed. 
This is because UGR with a cosmological constant admits an infinite number of 
maximally symmetric solutions labeled by the value of the space-time curvature. 
However, this does not explain why one should choose the desirable (almost) flat 
solution among a continuum of maximally symmetric ones. Another issue in UGR 



G{A)\V) = 



(1) 




(2) 



2 



is related to quantum loops of matter and gravity. It is likely that in this theory the 
Lagrange multiplier A acquires quadratic terms via the quantum loops (see, section 
3); if so, new interactions would be needed to maintain the classical solutions of 
UGR in a quantum theory. 

Constrained GR improves on the above-mentioned aspects: (i) it admits only two 
maximally symmetric solutions - one with zero curvature, and another one with the 
curvature obtained in GR; (ii) its classical properties, due to the BRST invariance 
of that model, are not modified by the quantum loop effects. 

Reparametrization invariance in CGR is completely constrained. We will show 
in section 3, that one can still define a locally inertial reference frame in a small 
region around an arbitrary space-time point. This is because the constraint that 
fixes gauge in the whole space, allows, in the neighborhood of a given point, for 
the point-dependent gauge transformations that locally eliminate effects of gravity. 
Thus the equivalence principle is preserved. 

Do the constrained theories solve the "old" cosmological constant problem [7j? 
In UGR the answer is negative because the theory admits an infinite number of 
maximally symmetric solutions. Clearly, in CGR, there still exists a conventional 
de Sitter solution of GR which can be used for inflation in the early universe. But 
the question is whether there exists an infinite number of other non-maximally 
symmetric solutions in the theory with a cosmological constant. If some of these 
solutions are physical, one should understand why in our Universe the (almost) 
flat solution is preferred. If, on the other hand, the non-maximally symmetric 
solutions can be disregarded for one reason or other, then CGR could be a good 
starting point for trying to accommodate inflation in the early universe and still 
solve the "old" cosmological constant problem. These important issues, including 
the question of stability of the new solutions, are being studied in jH]. The purpose 
of the present work is to investigate CGR to determine whether it is a consistent 
low-energy quantum field theory, which by itself is a legitimate theoretical question. 
Concrete applications of this model will be discussed in jS]. 

Even in the most optimistic scenario for CGR, one should explain why the space- 
time curvature is not exactly zero but ~ Hq ~ (10~ 42 GeV) 2 , as suggested by recent 
observations. Where could this scale come from? One possibility is to introduce a 
pseudo Nambu-Goldstone boson potential [0] that gives rise to required "dark en- 
ergy" 1 . Another way is to try introducing a graviton mass m g ~ H . Although 
both of the above approaches postulate the existence of a new small scale, this scale 
is stable w.r.t. quantum corrections (i.e., it is technically natural). In this regard, 
we briefly discuss in the present work massive deformations of classically constrained 
gauge and gravitational theories. We will find that the UV behavior of the prop- 
agators of massive gauge and gravitational quanta are softened. Although these 
results are encouraging, at this stage we still lack an understanding of whether the 
non-linear unitarity of the S-matrix on a Hilbert space of physical states can be pre- 

1 In this case though, one needs a VEV of a scalar to be somewhat higher than the Planck mass. 
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served using the BRST and ghost charges of these models, or whether some further 
modifications may be needed. Studies on these issues will be reported elsewhere. 

The work is organized as follows: in section 2 we discuss a constrained theory of 
a photon and its classical equations and solutions. In section 3 we discuss a classi- 
cally constrained theory of gravity (CGR). We find that CGR with a cosmological 
constant has a remarkable property - it admits two maximally symmetric solutions 
one of which is a flat space. We find general cosmological solutions in CGR as 
well as the expression for the Schwarzschild metric. An important question that is 
also addressed in section 3 is that of radiative stability. Using the BRST invariant 
version of the theory we argue that quantum corrections do not ruin the obtained 
classical results. Massive deformations of CGR are also briefly discussed in section 
3. In Appendix A we discuss the spectrum of the constrained Abelian gauge theory 
in various approaches, including Stiickelberg's method. We also look at the massive 
deformation of this model pointing out that the massive propagator, unlike in the 
Proca theory, has smooth UV behavior and a nonsingular massless limit. Appendix 
B deals with constrained non- Abelian gauge fields. After briefly discussing classical 
equations, we study the spectrum of the theory. Comments on massive deformations 
of the non- Abelian theories are also included. The results are similar to those of the 
Abelian case. Some parts of Appendix A and B are of a review character, but we 
felt that including these discussions would make our presentation more complete. 



2 Constrained photon 

As an instructive example we consider electrodynamics with an imposed classical 
constraint. We call this model constrained QED (CQED) even though we will only 
quantize it later. We start with the Lagrangian density 

/ = -^F + + \(d,An ■ (3) 

Here v — 0, 1, 2, 3; J M is a current, and A is a Lagrange multiplier (our choice of 
the Lorentzian signature is "mostly negative"). 

We start by studying classical properties of this theory. The equations of motion 
that follow from the above Lagrangian are 

& , F vlt + Jp-dp\ = 0, (4) 
V^ = 0. (5) 

Taking a derivative of Eq. (pEJ) one obtains the following relation 

□A = <9%. (6) 

If the current J M is conserved, A is a harmonic function, DA = 0. One particular 
solution of this equation is <9 M A = C M , where denotes an arbitrary space-time 
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constant four-vector. Physically, the effective conserved current to which the gauge 
field is coupled in J3J) is 



J'/f — Jfi~ d^X — Jp — Cfj,. (7) 

The last term on the r.h.s. of (J7|) acts as a constant background current density 
determined by a vector C M . Since C M is an integration constant, there is a contin- 
uous set of C M 's that one could choose from. Setting the value of is equivalent 
of choosing the corresponding boundary conditions. The equations of motion (J3J) 
and (0) were derived by varying the action corresponding to Q with the following 
boundary conditions 

(boundary = 0, 5X | boundary = finite function . (8) 

Therefore, all the solutions should obey (jSJ). To demonstrate that such solutions 
exist let us consider a simple example of a spherically symmetric localized charge 
density for which J = P#( r o — r), Jj = 0. Here p is a constant charge density, 
8(r) denotes the step function, r is the radial coordinate, and ro is the radius of the 
charge distribution. We substitute this source into the RHS of (jlj). In addition we 
chose a solution for the Lagrange multiplier to be d X = C = p, and djX = Cj = 0. 
For this source a solution of Eqs. (jlj) and (j5J), and the corresponding electric field 
read 

A) = - - for r > r , A = for r < r , 

Sr 6 2 

E = t (1 _ l|) 0{r - r )f . (9) 

The above solution satisfies the boundary conditions (JHJ at r = ro- In the conven- 
tional electrodynamics the source Jo = p0(r — r), Jj = 0, yields an electric field 
that is well-known and differs from (0). The origin of this difference is clear - in 
CQED the quantity J M that is specified in the action is not the whole source pro- 
ducing the gauge field! An additional integration constant appears in the equations 
of motion and the total source is J^ ff (|7jl. The above theory reduces to conventional 
electrodynamics when we choose A = 0. This corresponds to what we measure in 
ordinary experiments. 

Similar solutions with a nonzero value of the Lagrange multiplier will play an 
important role for gravity with a nonzero cosmo logical constant. These will be 
discussed in section 3 (non-Abelian gauge fields are discussed in Appendix B). 

So far we have not emphasized the fact that the Lagrangian Q is not gauge 
invariant. In fact, variation of Q under the gauge transformation 5A^ = d^a(x) 
vanishes (up to a surface term) for configurations satisfying DX = 0, therefore (J3J) 
has gauge invariance when the A field is on- shell, even though the gauge field could 
be off-shell. Given that the A is not propagating, and remains such at the quantum 
level (see below), this suggests two physical degrees of freedom for a photon, while 
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off-shell there are four degrees of freedom. This will be established more rigorously 
below. As it is shown in Appendix A, the two extra off-shell degree of freedom are 
decoupled from conserved sources and have no relevance in the Abelian case. 

The quantum theory of CQED could be approached in a number of different 
ways. One could restore first a manifest gauge invariance of Q using the Stiickelberg 
method, and then quantize the resulting gauge invariant theory. We will discuss this 
in Appendix A and B for Abelian and non- Abelian gauge fields respectively. 

In the Feynman integral formulation one could think of the constrained approach 
as follows: the gauge and auxiliary fields can be decomposed in their classical and 
quantum parts, A = A d + 5 A, and \ = \ c i + 5\. For classical solutions, A d and 
X c i, we allow boundary conditions that are different from the conventional ones. 
In particular, we allow for a nonzero solution of the DA = equation, a zero- 
mode. This is an unconventional step, since the usual FP term that appears in 
the path integral is the determinant of the operator acting on A, that is det{\3) 
in this particular case. Because of the zero-mode, the determinant, det(D), would 
have been zero, leading to an ill-defined partition function. However, a right way to 
formulate the path integral is to separate the zero-mode, and integrate only w.r.t. 
the small fluctuations for which det{\3) is non-zero, and for which the conventional 
radiation boundary conditions are imposed. 

This is the procedure that we will be assuming throughout the text. Then, 
quantization of the fluctuations 6 A and 5 A over the classical background A c i, A c j 
can be performed in a BRST invariant way. That is, one could start by postulating 
BRST invariance of the Lagrangian for the fluctuations by adding in the Faddeev- 
Popov (FP) ghosts without any reference to a local gauge symmetry, but elevating 
BRST invariance to a fundamental guiding principle in constructing the Lagrangian. 
The resulting BRST symmetric theory of small fluctuations could be quantized in a 
conventional way. This is what we summarize below. For simplicity of presentation 
we will replace 5A^ and 5X by A^ and A keeping in mind that these are fluctuations 
over a classical background (the same replacement will be assumed for non-Abelian 
and gravity fields considered in the next sections and Appendices). 

In a conventional Feynman integral approach the measure in the path integral 
should be modded out by the gauge equivalent classes. The FP trick does this job 
by introducing the gauge-fixing term along with the FP ghosts. In this regard the 
following natural question arises - since Q (or its non-Abelian counterpart) is not 
gauge invariant why do we need to introduce the FP ghosts in the theory? Naively, 
it would seem that we should use the path integral 




(10) 



instead of the one with the FP ghosts 




(11) 
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where G is a gauge- fixing condition (for instance it could be that G = d^A^), and 
G u is gauge transformed G with u being the transformation parameter (the above 
definition is consistent as far as we do not include the zero modes in the integration 
and det(n) is nonzero, as we discussed above). In QED the difference is irrelevant 
because the FP ghost are decoupled from the rest of the physics, however, in the 
case of non-Abelian fields (fTUJl would lead to a non-unitary theory. Because of the 
absence of gauge invariance in (j3J), the argument to introduce the FP ghosts in the 
path integral of (j2J) cannot be the same as in the conventional case. Nevertheless, 
the FP ghosts can be motivated by a symmetry argument, namely by requiring the 
BRST invariance of the constrained action. 

Let us modify Q by adding the FP ghost c and anti-ghost c which are Grassmann 
variables, c 2 = (c) 2 = and c + = c, c + = c. The Lagrangian reads as follows 

I = -\f^ + A^ + \(d»A») - icUc. (12) 

Since on the classical backgrounds considered here the FP ghost fields vanish, the 
classical properties of (JHJ) and (fT2*|) are identical 2 . The Lagrangian ()12j) . however, is 
invariant under the following continuous BRST transformations: 

5A^ = i(d^c, <5c = 0, 5c = AC, S\ = , (13) 

where ( is a coordinate independent Grassmann transformation parameter such that 
(C c ) + — C+ C + — C C- The Lagrangian (fT2~j) gives the path integral that is identical to 
that of the conventional approach (jllj) which takes the form 
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F, U F^ + A^J^ + \{d^) - icUc 



(14) 



Notice that the rationale for writing down ()14j) and the Lagrangian ()12j) is different 
from the motivation that led to the path integral (fTTj) . In the constrained approach 
the rules for constructing a Lagrangian and path integral are: 

(1) Impose classical constraints on fields using the Lagrange multiplier tech- 
nique. Find the corresponding zero-modes and treat them separately from small 
fluctuations of the fields. 

(2) Introduce the FP ghosts to obtain the BRST invariance of the gauge non- 
invariant theory. 

(3) Use this Lagrangian to set up the path integral in a straightforward way. 
The presence of the BRST symmetry guarantees that all the Ward-Takahashi 

identities (the Slavnov- Taylor identities in the non-Abelian case) of the conventional 
theory are preserved in the constrained approach, even though the classical equations 
of motion in this approach are different as discussed above. 

2 In general, a ghost-antighost bilinear could have nonzero expectation values on certain states, 
however, these stats do not satisfy the zero BRST and Ghost charge conditions, see below. 
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Let us now turn to the loop corrections that emerge in (|12|). In particular we 
would like to make sure that no kinetic or potential terms are generated for A. We 
do have a symmetry X(x) — > X(x) + f3(x), where (3 is an arbitrary function, w.r.t. 
which ()12j) is invariant. Kinetic or potential terms of A would break it. The question 
is whether this symmetry is preserved by the loop corrections. To address this issue 
we calculate the propagator of the gauge fields. This can be done in a few different 
ways. The easiest one is to add a fictitious term — \^X 2 to the Lagrangian density 
and then take the limit 7^0 



Z[JJ oclim / VAVXVcVcexp 
7-0 J 



oc lim / VAVcDc exp 
7—0 ' 



l g + A^J^ + X(d^) - ^ 7 A 2 



j lg + A% + ^(<V) 2 - ICUC 



icdc 



(15) 



where l g = -\F^F» V '. The propagator is obtained following the standard procedure 
and the result is 



A 



— lim 

7-+0 



jf - (1 + 1 )d^d v /U if - d»d v /U 



□ 



□ 



(16) 



which coincides with the standard transverse QED propagator in Landau gauge 3 . It 
is straightforward to check that loop corrections preserve the transversality of the 
gauge field propagator. All the diagrams that renormalize the propagator consist of 
the standard bubble diagrams produced by the contractions between two currents. 
Because the current to which the gauge field is coupled is conserved, the propagator 
remains transverse to all loops 



d»X (Ap A v )d v X ~ arx^u - d^d u )d v X = . 



(17) 



As a result, loops cannot generate the kinetic term for A. 

Last but not least, the BRST symmetry of (|T2|) can be used to define the Hilbert 
space of physical states by imposing the standard zero BRST- and Ghost- charge 
conditions Q BRST |Phys >= 0, Q Ghost |Phys >= 0. In the Abelian case this reduces 
|lj to the Nakanishi-Lautrup condition A^lPhys >= [H E] (here A 1- -* denotes 
a negative frequency part of a fluctuation of A over a classical background) which 
ensures that the quantum of A is not in a set of physical in and out states of the 
theory. 



3 Here and below we always assume the Feynman (causal) prescription for the poles in the 
propagators. 
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3 Constrained models of gravity 



There exists a modification of GR, unimodular gravity (UGR) j^j, that partially 
restricts its gauge freedom and yet reproduces the observable results of the Einstein 
theory. In UGR reparametrization invariance is restricted to the volume-preserving 
diffeomorphisms that keep the value of yjg intact (g = \detg^ v \). Such a theory 
can be formulated by using the Lagrange multiplier A (we set SnG^ = 1, unless 
indicated otherwise): 

C = -^R-\( y fg-l) + C M , (18) 

where Cm is the Lagrangian of the matter fields. The equations of motion and 
Bianchi identities of this theory require A to be an arbitrary space-time constant jH] • 
As a result, UGR is equivalent at the classical level to GR except that cosmo logical 
term becomes an arbitrary integration constant. Since the latter can take any value, 
there are an infinite number of solutions parametrized by a constant scalar curvature. 
Such a theory, at least at the classical level, seems to be more favorable than GR - 
for an arbitrary large value of the vacuum energy in the Lagrangian (arising, e.g., 
from particle physics) one is always able to find an almost-flat solution. Of course 
this still does not explain why one should choose the desirable almost-flat solution 
among a continuum of maximally symmetric ones. This is one aspect of UGR on 
which one would like to improve. 

Another, perhaps more pertinent question in UGR emerges when one considers 
quantum loops of matter and gravity. Inspection of loop diagrams (see section 4.2) 
suggest that the Lagrange multiplier A in (fTSjl would acquire the mass (A 2 ) as well as 
kinetic ((<9 M A) 2 ) terms due to the quantum effects. As a result, A would cease to be 
an auxiliary field, and all the classical results of UGR would have to be reconsidered. 

We will discuss a model that completely constrains reparametrization invariance 
of GR. This theory has the following two important properties: (I) It allows, like 
UGR does, an adjustment of the cosmological constant via the integration constant 
mechanism, but only admits two maximally symmetric solutions, one of which is a 
flat space. (II) Unlike UGR its classical properties are stable w.r.t. quantum correc- 
tions, i.e., the Lagrange multiplier remains an auxiliary field even in the quantum 
theory. This theory also preserves the equivalence principle. 



3.1 Constrained gravity and the cosmological constant 

Consider the following Lagrangian 

C = ~^yR + y/gg^d^\ v + C M + surface terms. (19) 

Here, A„ is a vector that serves as a Lagrange multiplier and Cm denotes the La- 
grangian of other fields which can also include a vacuum energy term (the cosmo- 
logical constant) produced by classical and/or quantum effects. Versions of this 
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model have been discussed in the literature previously (see, e.g., the last reference 
in [3]) with the purpose of introducing the de Donder gauge- fixing condition in the 
context of quantization of gravity. Here, instead, we regard this model as a classical 
theory, which is subsequently quantized, but the classical equations of which could 
admit rather interesting solutions that are absent in GR. Thus, we first concen- 
trate on classical effects, leaving the discussion of quantum corrections for the next 
subsection 4 . 

Among others, we will be considering below solutions with fixed boundaries. For 
such solutions one should add to the Einstein-Hilbert action the Gibbons-Hawking 
boundary term. Moreover, the boundary conditions that we allow for are: 

^9imu I boundary = 0, and SX^ | boundary = . (20) 

Under these conditions the variation w.r.t. the metric gives 

Giw - (<9A + d »K) + 9^K = . (21) 

In addition to this, we get a constraining equation by varying the action w.r.t. the 
Lagrange multiplier 

^{y/g g^) = ■ (22) 

The above model is similar in spirit to a constrained gauge theory discussed in the 
previous section. It is easily checked that (f22~j) is equivalent to the condition 

v^r = 0, (23) 

and, in the linearized approximation, gives rise to de Donder (harmonic) gauge fixing 
of linearized GR. Due to this condition, V^A^ = g^ u d^A u = d^A^, where is a 
covariant derivative acting on an arbitrary four- vector A^. 

The constraint (J22|) (or (J23j0 fixes completely reparametrization invariance of the 
theory. How is then the equivalence principle recovered? A relevant property of (J22|) 
(or (JHU)) is this: for any given point x^ = Xq in the coordinate system {x 11 }, it allows 
for the ^Q-dependent coordinate transformations that eliminate the connection in a 
small neighborhood of this point. These transformations can be written as 

z" 1 = x» + l -T^{x ){x - x ) a (x - x,f . (24) 

It is straightforward to check that for (J2"4~j) g a/3 (x)\ x=X0 = g' al3 (x')\ x i =X0 , and 

^^{ X )\x=x — r^(x )\x'=x + ^a/3( X )\o==x ■ (25) 

4 One could choose to impose a different constraint using the Lagrange multiplier in (|f 9fl (for 
instance, an axial-gauge constraint). As long as it is an acceptable gauge-fixing condition for small 
fluctuations, and the corresponding FP ghosts are taken care of consistently (in the axial-gauge 
the FP ghost are not needed) the quantum effects of the fluctuations won't depend on the choice 
of this constraint. However, the classical solutions could differ for different constraints. Here we 
choose the constraint that is Lorentz- invariant. 
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As a result, r^ 3 (x / )| x / =Xo = 0, and at this point the metric can simultaneously be 
brought to the Minkowski form. The transformation (|24|) . on the surface of (J22)) (or 
dZflO, is trivially consistent with (I22J) for lj25 |l ) 

^I^COU* = ^(x')O(x')|^=, = 0. (26) 

Since x was arbitrary, the above arguments can be repeated for any other point in 
space-time. One should emphasize again that the coordinate transformation (j2H 
is point-specific, i.e., at different points of space-time one should perform transfor- 
mations that depend parametrically on that very point. This is why they are not 
gauge transformations in the entire space-time. Summarizing, although (J22j) (or 
(j22J)) picks a global coordinate frame, it allows for the point-dependent coordinate 
transformations that can eliminate gravity locally. Hence, the equivalence principle. 
The local Lorentz transformations are also preserved. 

To obtain the equation which the Lagrange multiplier has to satisfy we apply 
a covariant derivative to both sides of (f2"Tj) . This is subtle since the second term 
on the l.h.s. of (|2"T]) is not a tensor and we should define the action of a covariant 
derivative on this object. We adopt the following straightforward procedure: apply 
to both sides of (|2"Tj) the operator 

V a = gT (d, - - r:j , (27) 

where the standard index arrangement should be used in place of the asterisks even if 
the two-index object on which this operator is acting does not transform as a tensor. 
Then, using the Bianchi identities and covariant conservation of the stress-tensor we 
obtain: 

g^^d a X u + d v K) = V^ CT A ff . (28) 

The above equation can be simplified substantially due to (f2~2"j) . The left hand side 
of pSj) can be reduced to 

g m d» {d a X v + d v X a ) - (d a X + d p X a ) , 

while the right hand side simplifies to give 

{d v g aP )d a Xp + g aP d v d a Xp. 

Combining the above two expressions together we find from (|2"8*j) 

g^ a 9fj,d a X v = . (29) 

As long as g^ v is non-singular general solutions for A„ could be found. 

Clearly, the system of equations (|2"Tj) . (}2"2"|) and (J2"8~j) (or (}2"9*j0 . could admit new 
solutions that are absent in the Einstein theory. For instance, the Einstein equations 
with a nonzero cosmological constant (i.e., T^ v = Ag^ u ), do not admit Minkowski 
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space as a solution. In contrast with this, the system (|21|). ()22|) and ()28|) is satisfied 
by the flat space metric 



S> = c^A,, + c^A^ = A 77^ . (30) 

It is remarkable that one can obtain a flat solution even though the vacuum energy in 
the Lagrangian is not zero! Similar property exists in unimodular gravity [U] when 
T^y = Ag^, but there are an infinite number of maximally symmetric solutions 
labeled by the value of a constant curvature. Is this also true in the model (fT^j)? We 
analyze this issue below. First we notice that the system (}2*Tj) . (j2^j) and (J2~%j) admits 
a de Sitter solution of conventional general relativity (we will focus on the case of 
a positive cosmological constant only from now on and interesting results can be 
obtained for both positive and negative A, jH]). The conventional dS metric solves 
(|2*T|) with A^ = const.. What is less obvious is how this solution satisfies (j22j) . To 
understand this we start with the dS solution in the co-moving coordinate system 
(this is also applicable to any FRW cosmology) 

ds 2 = dt 2 -a 2 (t)(dx 2 + dy 2 + dz 2 ). (31) 

As it can be checked directly, (|3*T| does not satisfy ((221) ■ However, we can define a 
new time variable 

f dt 

T -IWY (32) 

for which the interval becomes 

ds 2 = a\r)dr 2 - a 2 (T)(dx 2 + dy 2 + dz 2 ). (33) 

This metric satisfies (|22j) identically. Therefore, a conventional dS metric, or any 
other spatially-flat FRW cosmology, is a solution of the system (12 1|) . (j22J), (EH)- 

Having established that the flat and conventional dS spaces are two solutions 
of the theory, let us now look at other possible maximally symmetric solutions. In 
general, the following ansatz 

9A + ^ M = <, (34) 

where c is an arbitrary constant, gj£ is a dS solution with R = — 4(A — c), does 
satisfy the equations (|2Tj). (J22|) and (|25|) . However, this is not enough to claim that 
(J53)l is a legitimate solution of the theory. This is because equation (jMj) itself might 
not be solvable in terms of A^ given that is a dS metric obeying (J22J); solvability 
for A M is a necessary condition since it is w.r.t. A M that we varied the action. It is 
straightforward to check that there is no solution for A M that would satisfy equation 
(jMj) if the metric is given by Could there be other forms of dS space that are 
non-trivially different from (|33|) and yet satisfy (|22j) ? The answer is no. To see this 
consider a dS solution in GR in two different coordinate systems. Let us assume the 
opposite, that both of these coordinate systems can be gauge transformed in GR to 
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two different coordinate systems for which fl2H) is valid. If true, this would mean 
that the condition (J2~2")) does not completely fix the gauge freedom of GR. On the 
other hand, we find by performing gauge transformation of (|22j) that this is only 
possible if the gauge transformation itself is trivial. Therefore, the form (J3*3j) is the 
unique dS solution which satisfies the constraints (J221)- 

The fact that (j34j) cannot be a maximally symmetric solution with nonzero R(g) 
could also be established just by looking at a general expression of the Ricci scalar 
in terms of the metric g^ v in which the substitution cg^ v = d^X u + d u \^ is made. If 
this ansatz could describe a dS space, one could always go to a weak-field regime 
where the Ricci curvature, as a function of the metric, has to be nonzero. How- 
ever, the above ansatz for g^ v gives zero R in the leading order of the weak- field 
approximation. 

Summarizing, we conclude that in the class of maximally symmetric spaces there 
are only two solutions of the theory: (I) The flat space defined in (j3l]j) ; (II) The (A)dS 
solution as it would appear in conventional GR transformed to a new coordinate 
system. There could in principle exist other, non-maximally symmetric solutions. 
It would be interesting to study whether those solutions are physical. If they are 
one should look for arguments why the maximally symmetric solutions could be 
preferred in our Universe. On the other hand, if the non-maximally-symmetric 
solutions are not there, then the model (fT^j) could be a playground for studying 
the fate of the "old" cosmological constant problem |7]. It is interesting to point 
out that by integrating out the Lagrange multiplier field one gets a non-local action, 
with some similarities but also differences of the resulting non-local theory with that 
ofRef. HD]. 

We would also point out that the Schwarzschild metric of conventional GR is 
also a legitimate solution of the CGR. A simplest way to address this is to choose 
X u = and make sure that the known GR solution itself satisfies (J2~2"|) in a particular 
coordinate system. Usual solutions of GR can be transformed to satisfy Eq. (f22|) . 
This is in particular true when g^ v is diagonal and each element of ^fgg^ (there 
is no summation w.r.t. ji here) factorizes into the products of the form ^fgg^ = 
h(x a )j(x u ) ■ ■ ■ f(x x ), where each function depends on one coordinate only. In this 
case, the constraints ()22|) turn into four separate partial differential equations 



Suppose we introduce one new coordinate / = x a (x a ) such that it depends only 
on x a , and leave all the other coordinates intact. In the new coordinate system 



/i = 0, 1, 2, 3, no summation w.r.t. fi. 



(35) 




(36) 



while 



w , and g^ = g^ for fi ^ a, 
and chosen it to be positive. 



(37) 



where we have defined x' a = 
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The a-th equation of ()35|) in the new coordinate system takes the form: 



dx a dx a ' ^ ' 

If y/—gg aa can be factorized, say as, \/—gg aa = h(x a )i/j where ip depends only on 
coordinates other than x a , we can find the desired x a by simply demanding 

h(x a ) 

-— ■ = 1 , or any constant if more convenient, (39) 



and solving this ordinary differential equation. It is not difficult to see that one can 
carry on the same procedure for each x^'s without invalidating the previous results, 
and, therefore, eventually find the new coordinate system that satisfies (}2*2*|) . 

The above procedure is directly applicable to the Schwarzschild and FRW solu- 
tions. For the latter the result was already given above (see (|33Jl ). Here we perform 
the change of coordinates for the Schwarzschild metric. In a spherically symmetric 
coordinates 

ds 2 = (l - -) dt 2 - (l - -) 1 dr 2 - r 2 (d6 2 + sin 2 6d<p 2 ) . (40) 

Here, r g = 2G^M is the horizon radius of an object. The above described procedure 
leads to the new coordinate system 



r - r g r g 
r g m =>- r — 



1 - e f / r 



3 



9 

6 = In tan- = 2tan~V. 



(41) 



2 

in which the Schwarzschild metric becomes 



ds 2 = e f ^dt 2 - —^——-dr 2 T i ( d6 2 + dip 2 ) . (42) 

(e r ^-l) 4 { e rlr 3 - 1)2 cosh 2 Q \ ^ > 

This metric satisfies ()22|) . The new variable f spans the interval (— oo, 0] as the 
coordinate r increases from r g to +oo, while the new angular variable 9 covers the 
interval (— oo, + oo). One can also easily describe the interior of the Schwarzschild 
solution by flipping the sign of the argument of the log in ()41|1. 



3.2 Radiative stability 

In this section we discuss the issue of quantum loop corrections to ([T9*|). In particular, 
we would like to understand whether this theory is stable w.r.t. the loops. It is 
clear that quantum-gravitational and matter loops will generate higher dimensional 
operators entering the action suppressed by the UV cutoff of this theory. This 
is similar to any theory that is not renormalizable and should be regarded as an 
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effective field theory below its UV cutoff (for an introduction to an effective field 
theory treatment of gravity, see, e.g., PI]). However, there is another question of 
a vital importance in the present context. This is whether loop corrections can 
generate the potential and/or kinetic terms for the Lagrange multiplier. If this 
happens, A M cannot be regarded as an auxiliary field and all the results of the 
previous subsection would be ruined. 

We will argue below that this problem is avoided in CGR because of the specific 
form of (|19j) which can be completed to a BRST invariant theory. Following Ref. 
[T2] , we introduce the variables (as in the previous sections, below we are discussing 
small fluctuations of the fields) 

-r = y/g<r, 7^ = ^- (43) 

It is straightforward to rewrite the Lagrangian (|T9*|) in terms of these variables and 
include the FP ghost term: 

C = - l -Y u 0M7) " d,K - 3,A M ) + l - (d,c u + 0„c„) Vf c a . (44) 

Here the terms in the first parenthesis represent the gravitational part of ([19)1 . while 
the last term introduces the vector-like FP ghost and anti-ghost fields for which the 
operator is defined as follows: 

The important point is that the standard Einstein-Hilbert action can be rewritten 
in terms of (|43jl and their first derivatives [T2j 

- \^ V R,M = -\d pl ^d al x » (i^^rv - \l pa l,rlxu - 2«7m*) • (45) 

The FP ghost term in (|44|) ensures the BRST invariance of this theory ^3]- The 
respective BRST transformations with a continuous Grassmann variable ( are: 

5^ = ic T d T c»(, 
5c» = -A M C, <^V = . 

Presence of this symmetry ensures that the Lagrange multiplier in (}4*4"|) does not 
acquire the kinetic term through the loop corrections. This is because the potentially 
dangerous term in the effective Lagrangian 

^A^Orh^O))^, (46) 

is zero (up to a total derivatives) due to the transversality of a two point graviton 
correlation function. The latter can be seen order by order in perturbation theory. 
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Figure 1: The solid lines denote the A field (there is no propagator corresponding to these 
lines, they are depicted to show the vertices); the wave lines correspond to gravitons. The 
vertices in this one-loop diagram arise due to the cubic interaction of A with two gravitons 
originating from the term Xy/g in the Lagrangian (j!8|) . 

Let us for simplicity consider this for an expansion about a flat space which is a 
consistent solution of the theory even if the cosmo logical constant is present. We 
introduce the notations 7 Miy = r]' iu — and look at the two-point correlation func- 
tion of the ip^ v field. This correlator has been studied in detail in the conventional 
approach, in which there is no Lagrange multiplier term in (|44j) . but instead, a stan- 
dard quadratic gauge fixing term ) with the gauge parameter ( is introduced. 
Our theory, on a fixed classical background, corresponds to the limit ( — > 0. Hence, 
all the results concerning the quantum loops derived in the conventional approach 
on a fixed background are also applicable here with the condition that ( — > 0. The 
BRST invariance of the theory can be used to deduce the Slavnov- Taylor identities 
[T4"] (see, e.g., [THJ E3 El)- The latter guarantee that order by order in perturba- 
tion theory the two point correlation function of the ^ v field is transverse (this 
corresponds to the ( — > gauge results of Refs. [13 CHI El) 

(<p^(x)<p al3 (o)) oc n^n 1 ^ + u^u ua - iptt^, 

where IP" = rf v . 

In conclusion, the presence or absence of the FP ghosts does not affect the classical 
equations of motion, and for the classical analysis it is acceptable to ignore them 
and study the Lagrangian ()19p . The only quantum mechanically consistent theory 
is that with the FP ghosts described by ([4"4"j) and all the classical results obtained 
above hold in this theory. Furthermore those results are stable w.r.t. the quantum 
loop corrections. 

We would like to comment on a similar issue in the context of UGR [§]. The 
Lagrangian of this theory, as given in ()18|). is likely to generate quadratic terms for 
A via loops. To see this we start with a one- loop diagram of Fig. 1. This diagram 
is logarithmically divergent and will generate an additional term proportional to 

^log {^A , (47) 

M P1 V »IR J 

where \i\jy and [im denote the UV and IR scales respectively, and we have restored 
Mpj in front of the which resulted in the 1/Mpj coefficient in ()47|) . Naively, the 
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Figure 2: The solid lines denote the A field (there is no propagator corresponding to these 
lines, they are depicted to show the vertices); the wave lines correspond to gravitons. 
The vertices in this one-loop diagram arise due to the quartic interaction of A with three 
gravitons originating from the term X^/g in the Lagrangian ()18|) . 

above term may seem to be irrelevant because of the 1/Mpj suppression. However, 
to get a right scaling we should restore the canonical dimensionality of A. This 
is achieved by substitution A — > Mp^a, after which the a-dependent terms in the 
effective Lagrangian take the form 

M&iy/g - 1) + M^a 2 log . (48) 

Because of the new induced term, the field a acquires a Planckian mass and ceases 
to be an auxiliary field. Furthermore, we could also look at a two-loop diagram of 
Fig 2. A simple power-counting of the momenta running in the loops shows that 
this diagram generates not only a mass term for a but also its kinetic term: 

Higher loops are also expected to generate similar terms, and the above arguments 
are hard to avoid unless the theory (JTSJl is amended by new interactions. Perhaps the 
BRST invariant completion of UGR proposed in Ref. can cure this; it would be 
interesting to perform explicit calculations in the framework of [17 to see whether 
the radiative stability is restored. Since the model of ^7j is BRST invariant one 
would expect a positive outcome. However, we should point out that the bosonic 
part of 17J contains additional gauge fixing terms needed to completely restrict 
parametrization invariance of the theory, and, from this perspective, it differs from 
UGR. 

Finally, we would like to calculate a response of the graviton field to a source. In 
the linearized approximation the number of physical propagating degrees of freedom 
of CGR should be the same as in GR. In the linearized approximation = 

due to (|2*2*|) . and d 2 X u = due to (|2*5jl . Then, the equation (}2~T]) simplifies to 

= . (49) 

The above equation is identical to an expression for the response in the Einstein 
theory. This of course is a consequence of the fact that the free propagator of 
(J44)) . coincides with the graviton propagator of GR in the harmonic gauge d ll ^ v = 
-dphT + \d v h = 0, where ~ if - h? v . 
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3.3 Comments on massive theories 



It is difficult to construct a consistent Lorentz-invariant nonlinear model of a massive 
graviton propagating on a Minkowski background. In the linearized approximation 
the only consistent massive deformation of GR is the Fierz-Pauli (F-P) model [TH| 

C = -^-R + ^m 2 g (h^ - h 2 ) , (50) 

were = g^ — rj^. In the quadratic approximation this model describes a massive 
spin-2 state with 5 degrees of freedom. However, a nonlinear completion of this the- 
ory is not unique, and so far there is no known non-linear theory in four-dimensions 
that would be consistent. A rather general class of nonlinear completions of the F-P 
theory give rise to unbounded from below Hamiltonian [T^]. This manifests itself in 
classical instabilities for which the time scale can be substantially shorter than the 
scale of the inverse graviton mass EQ 122] • 

The reason for this instability is that at the nonlinear level a ghost-like sixth 
"degree of freedom" shows up. This should have been expected because of the fol- 
lowing. Ten degrees of freedom of g^ v in (|5Uj) are restricted only by four independent 
Bianchi identities. Hence, six degrees of freedom should remain. The absence of the 
sixth degree of freedom in the linearized theory was just an artifact of the linearized 
approximation itself jjlGL 

It is interesting to ask the following question: what if we start with a mass 
deformation of a constrained gravitational theory instead of modifying GR as in 
(150(1 ? Straightforward calculations show that the mass deformation of unimodular 
gravity (|18j) or the constrained gravity (|19() leads to a theory with a ghost in the 
linearized approximation. This ghost can be removed, at least in the linearized 
theory, if one considers a mass deformation of a model with both constraints ^fg = 
1 and d^^fgg^ = 0. Because this set of equations imposes 5 conditions on ten 
components of g^ u , one should expect this theory to propagate 5 physical degrees of 
freedom. Below we will discuss the advantages as well as difficulties of this approach. 

The Lagrangian of the massive "hybrid model" that combines the two constraints 
mentioned above takes the form: 

£ = -^Lr+ y^ m l(h^ - h 2 ) + y/ggTdpK ~ HV9 ~ 1) • (51) 
Variation of the action w.r.t.the Lagrange multipliers A and A M yields the constraints: 

V£=i; <9 MV /^r = o. (52) 

We now turn to the linearized approximation about a flat space to study a 
graviton propagator. In this approximation the constraints (j52| reduce to h = 0, 
and d^h^y = respectively. The equation of motion becomes 

G,„, - m 2 g h^ - dfr\ v ) + {d a X a )r] IMU - Xr]^ = T^. (53) 
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The trace equation and the Bianchi identities give respectively 



-4A + 20% = T , d M A + <9% = , 
which can be solved to obtain 



d T T 
A» = ^ft, A = --. (54) 
6D 6 v ; 

Substituting these solutions into we find 

O + m* + 3D(D + m2)' 1 j 

From the tensorial structure of (jo3j) we conclude that the theory propagates five 
physical polarizations, as it should. Moreover, unlike F-P gravity, the propagator 
(J55)) has a well-defined m g —>■ limit. This is similar to the soft behavior of mas- 
sive gauge field propagator discussed in Appendix A and B, and to "softly massive 
gravity" emerging in higher dimensional constructions [2HJ E3] • 

Could (|51|) be a consistent model of a massive graviton in 4D? There still is a 
long way to go in order to find out whether (joTj) is a theoretically sound theory. 
There are three major checks one should perform. 

(i) The main problem of the F-P gravity stems from the fact that the Hamiltonian 
of the nonlinear theory is unbounded below. Hence, one should understand whether 
the same problem is evaded by the hybrid model (JoTj) . We studied this question 
partially and have shown that the unbounded terms that appear in the F-P massive 
gravity do not arise in (|51j). To understand this we look at the ADM decomposition 
of the metric 



In terms of the new variables the Einstein-Hilbert Lagrangian takes the form 

-^-N{R® + KijRV - K 2 ) , (56) 

where 7 = det^-, is the 3-dimension Ricci curvature calculated with the metric 
7ij and the extrinsic curvature tensor 

contains a covariant derivative Di with respect to the metric 7^. The problem of 
the F-P gravity arises because the lapse iV acquires a quadratic term in non-linear 
realizations of F-P gravity. Hence, it ceases to be a Lagrange multiplier and does not 
restrict the propagation of an extra sixth "degree of freedom" which is ghost-like. 
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In the hybrid model, however, because of Jg = 1 we get that N = 1/ >J\det ; y fll/ \. 
This constraint enables one to remove the sixth degree of freedom and the dangerous 
terms that previously led to unboundedness of the F-P Hamiltonian [T^j. Regretfully, 
the expression for the Hamiltonian of the hybrid model is rather complicated and it 
is difficult to see that there are no other sources of rapid classical instability there. 

(ii) Even if the rapid classical instabilities are removed, the main question is 
whether the BRST invariant completion of (J51|) exists. The BRST symmetry, could 
be a guiding principle determining a unique nonlinear completion of (j51j) for any 
other massive theory), which at this stage is completely arbitrary. One could hope 
that the BRST and ghost charges can be used to define the Hilbert space of physical 
states of the theory so that even if the Hamiltonian is not bounded below, the states 
of negative energy do not appear in the final states (i.e., they are projected out 
by the conditions QBRST|Phys) = QchostlPhys) = and cannot be emitted in any 
process). At the moment it is not clear whether such a construction is possible, but 
we plan to return to this set of questions in future. 

(iii) The question of radiative stability of (joTj) is something one should worry 
about. In general the Lagrange multipliers of the massive theory will acquire the 
mass and kinetic terms, and this would lead to propagation of a new degree of 
freedom. Only hope here could be to complete (J51|) in a BRST invariant way so 
that the resulting theory does not generate the quadratic and higher terms for the 
Lagrange multipliers. 

Even if all the above three issues (i— iii) are positively resolved, one needs to amend 
the massive model to make it consistent with the data. The point is that the scalar 
polarization of a massive graviton couples to sources and gives rise to contradictions 
with the Solar system data. In the model of Ref. the similar problems is solved 
due to nonlinear effects that screen the undesirable scalar polarization at observable 
distances [211 123 I2S EH] ■ In the present model, at least naively, such a mechanism 
does not seem to be operative, and some new ideas are needed. 
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Appendix A: Qantization of CQED: Stiickelberg formalism 



There is another way of quantizing We can restore the gauge invariance 
of Q using the Stiickelberg method and then follow the standard FP procedure of 
fixing the gauge and introducing the FP ghosts. Let us discuss this in some detail. 

We start by rewriting (jSJ) as follows: 

£ = -\rfr + A(0"B„ - Op) + J"(fl„ - d,<p), (57) 

here we have introduced the notations: 

B fl = A IJi + d^ip and = d^B u - d v B^ = d^A u - d u A^. (58) 

A Lagrangian, similar to ([57]). in a context of a theory with a non-conserved current 
was recently discussed in [3U]. The Lagrangian (|57|) is invariant under the following 
gauge transformation 

SB^ = d^a and Sip = a, (59) 

where a is an arbitrary function. To fix this freedom we choose a gauge similar to 
the i?£-gauge that eliminates a mixing terms between A and B. This can be achieved 
by adding into the Lagrangian the following gauge fixing term: 

A£ GF = ±(d»B,-Z\) 2 . (60) 

Here, £ is an arbitrary gauge parameter. Furthermore, it is easy to find the FP 
determinant and introduce the FP ghosts into the theory. The total Lagrangian 
that includes the gauge fixing and FP ghost terms reads: 

Act = -\f%, + hd^B^f - icUc - XOp + |A 2 + J»(B, - d^cp). (61) 

The first three terms on the right hand side of the above expression constitute a 
free Lagrangian of gauge-fixed QED. In addition there are other states in (|61jl . To 
understand their nature, we integrate out from (|61j) the auxiliary field A. As a result, 
the terms of (jfiTj) containing A get replaced as 

-AD^ + |a 2 -,-1(^) 2 . (62) 

Because the propagator of ip is ^ = lim^o^i^ — it is more appropriate 

to think of two states, described by (JH2J) one of which has a positive-sign kinetic 
term and the other one is ghost-like 5 . However, these field are not present in the 

5 This could also be understood directly from (|61J) by making substitutions: A = a + b, ip = a—b, 
which generate two kinetic terms one for a with a right sign and another one for b with a wrong 
sign. These two fields will also have ^-dependent masses and mass mixing that originate from the 
£A 2 term. 
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physical on-shell spectrum of the theory. In the limit £ — > the field tp is frozen, the 
Lorentz-gauge fixing condition, d^B^ = 0, is enforced, and the on-shell spectrum 
consists of two physical polarizations of a photon. Since physics cannot depend on a 
choice of the value of the gauge parameter £, the model (JoTj) propagates two on-shell 
polarizations. Off-shell, however, there are four propagating degrees of freedom. 
The two non-physical degrees of freedom are longitudinal and time-like components 
of a photon, while A plays a role of the canonically conjugate momentum to A 
(alternatively, if A is adopted as a canonical coordinate A becomes its conjugate 
momentum) . 

It is straightforward to read off the propagators from 

A g = _V li+jjg o /u ^ A ^ = _i^_. (63) 

Adding these two contributions together we find that the dependence on the gauge 
parameter £ cancels out and we obtain 1)16)1. This confirms our previous conclusion 
that the (p field is a gauge artifact. It is also instructive to rewrite the Lagrangian 
(JoT]) in terms of the original field A^. The latter looks as follows: 

Act = -\f% + ^A,) 2 - hPAJOp - icUc + A^ , (64) 

where we have integrated out the A field. Here we see a difference from conventional 
QED. The additional £ dependent term ensures that the resulting free propagator, 
for arbitrary values of £, coincides with the Landau gauge propagator of QED. This 
is due to an extra field ip which has no kinetic term but acquires one through the 
kinetic mixing with the gauge field 1)64)1 . The mixing term itself is gauge-parameter 
dependent and this is what cancels the £ dependence of the QED propagator. 



A.l. On massive deformation of CQED 

In this subsection we study the massive deformation of the Lagrangian Q. We 
expect, because of the classical constraint, the massive theory to be different off- 
shell from the conventional massive electrodynamics (the Proca theory). 
To proceed, we add to the Lagrangian (J3J) the following mass term: 

A/ = l -m 2 Al . (65) 

One can think of this term as arising from some higher-dimensional operator in which 
certain fields acquiring VEV's generate (J65)) while these fields themselves become 
heavy and decouple from the low-energy theory. 

It is straightforward to see that Eq. (jlj) gets modified by the mass term 

d v F vll + rrijA^ + J M - O^X = , (66) 
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while the constraint d^A^ = remains intact. Taking a derivative of ()66|) one 
obtains the equation of motion for A 



□ A = <9 M J M . (67) 

From the above we find a solution with a nonzero constant background current 
d^X = which is identical to that of the constrained massless theory. 

It is interesting that the constrained massive theory also has a continuous BRST 
invariance if the FP ghost fields are introduced. Indeed, consider the Lagrangian 

/ = -\f^ + \{d^) + ~m*Aj - icUc. (68) 



As before, the presence of the FP ghost in ()68|) does not affect the discussions of 
the classical equations. On the other hand, due to these fields the Lagrangian (jSHJ) 
becomes invariant under the following continuous BRST transformations: 

5A fl = i(d fi c, <5c =0, 5c=X(, 5X = im 2 (c, (69) 

where ( is a continuous Grassmann parameter. Notice that A transforms to com- 
pensate for the non-invariance of the mass term 6 . Moreover, 5 2 A fl = 5 2 c = 5 2 X = 0, 
while 5 2 c ^ 0, but 5 3 c = 0. 

Let us now turn to the discussion of the spectrum of this theory. First we evaluate 
the propagator of the gauge field. For this we follow the method used in section 2. 
The result is: 



A"" = - lim 

7^0 



1+7 
J_fQI*Q» 



□ — 7m, 



1 _ o»d u /n 



□ + m 2 □ + m 2 



(70) 



This should be contrasted with the propagator of conventional massive QED (the 
Proca theory): 

rT + d^/m 2 

A P-a- n + m 2 ■ {fL) 

The key difference of ()7())1 from (fTT|) is the absence in (fTT)|) of the longitudinal term 
that is inversely proportional to the mass square. Because of this, the propagator 
(J70J) does have a good UV behavior, while the propagator (fTTj) does not. In an 
Abelian theory with a conserved current this hardly matters since the longitudinal 
parts of the propagators do not contribute to physical amplitudes. However, this 
could become important for non- Abelian and gravitational theories where the matter 
currents are only covariantly conserved, so we carry on with this discussion. 

6 The above transformations H69(l differ from the standard BRST transformations of a massive 
Abelian theory in which massive FP ghost need to be introduced. We emphasize here the transfor- 
mations (|69fl because they will be straightforwardly generalized to the non- Abelian massive case 
in the next section. 
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A natural question that arises is what is the mechanism that softens the UV be- 
havior of (J7UJ) as compared to (|71|)? Usually in massive gauge theories this is achieved 
by introducing a Higgs field that regulates the UV behavior of (|71|). Therefore, on 
top of the three physical polarizations of a massive gauge field, there should be a 
new state that replaces the role of the Higgs. To see this state manifestly we rewrite 
(fTUj) as follows: 

□ + m 2 rrgp' 1 ' 

The first term on the r.h.s. is just the Proca propagator of three massive polariza- 
tions; the additional term represents a new massless derivatively-coupled degree of 
freedom. To uncover the nature of this extra state we rewrite the Lagrangian of the 
constrained massive theory as follows: 

l = + i™? {K " ^) 2 " " + W ■ (*>) 

Notice that the quantity in the parenthesis is invariant under the BRST transforma- 
tions (|69|h while the non-invariance of the A- kinetic term under ()69|) is compensated 
by the terms coming from the FP ghost kinetic term. Defining a field 

B» = A»- d -4, (74) 

we end up with the following theory: 

I = -\F% + \m*Bl + - ^(d»\) 2 - ^-c^c + (75) 

The first three terms on the r.h.s. of the above expression represent the Proca 
Lagrangian of massive electrodynamics. There are also additional terms in (|75|). 
These are the A kinetic term and FP ghost kinetic term. These two terms form a 
sector of the theory that is invariant under the continuous BRST transformations 
5X = im^(c, 5c = (X and 5c = 0. This symmetry is an exact one if the current 
is either conserved, as it is in an Abelian theory, or transforms w.r.t. BRST in an 
appropriate way as it will for a non-Abelian theory discussed in the next section. 
The A kinetic term in (|75j) has a wrong sign, and this state is ghost-like. While J M is 
conserved, A is decoupled from the rest of the physics and can be ignored for all the 
practical purposes. Nevertheless, it is interesting to understand whether the state 
A could belong to a Hilbert space of physical states. This space could be defined by 
introducing the BRST and Ghost charges and postulating Q BRST |Phys >= and 
<5 Ghost |Phys >= 0. The above conditions could ensure that quanta of the A field 
do not belong to the set of physical in or out states of the theory. However, as 
we discussed previously, the BRST transformations (|69|) are rather peculiar and the 
construction of the Hilbert space of physical states could differ from the conventional 
one. 
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Here we discuss the properties of the physical degrees of freedom in the Hamil- 
tonian formalism. Let us ignore the external current J M for the time being. The 
canonical momenta conjugate to read 

dl 

The key difference of this from QED is that (j2J) contains time derivative of A , and, 
as a result, the primary constraint of QED, tt = 0, is replaced by the relation 

tto = A. (77) 

Thus, A is just a conjugate momentum for A . The above expression can be used to 
determine A, and if so, it does not constrain n . Excluding A by (|77jl . the extended 
Hamiltonian takes the form 

H = -ttJ + ^(tijkdjAk) 2 + TTodjAj + A (djTTj) . (78) 

To study this system further we look at the two equations governing the time evo- 
lutions of 7r and A as 

ir = {n H}, A = {A H}, (79) 

where {..} denotes the canonical Poisson brackets. This reduces to the following 
relations 

7T + d^i = 0, d^A^O, (80) 

both of which were already given in a different form by equations of motion derived in 
section 2. Furthermore, requiring that the time derivatives of these two expressions 
are identically zero 

{7r + d&i,H} = 0, {d^,H} = 0, (81) 

we find the following two additional equations 

□ 7t = , □ A° + diTTi = 0. (82) 

Further time derivatives are identically satisfied. 

From the first equation of (|K2*|) we find that ttq can either be a plane wave or a 
trivial harmonic function (or a superposition of the two). The second equation of 
([82)1 . however, dictates that ttq can only assume the trivial solutions. To show this 
suppose that the solution for n is a plane wave. Then, according to 

DA = -dm = 7T , (83) 
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A cannot have a nonsingular solution because of the mass-shell condition. There- 
fore, to avoid non-physical solutions a trivial solution for ir should be taken. As a 
result, equation (|8U|) and ([82]) together remove 7r and A from the list of free vari- 
ables and impose one condition on tt^s and one on A l 's, thus, reducing the number 
of physical degrees of freedom to two. This is consistent with the Nakanishi-Lautrup 
condition that prescribes to physical states to satisfy A^|Phys) = ttq ^|Phys) = 0. 

In the massive case the Hamiltonian density is modified to be 

n = \v) + \{e l3k d 3 A k f + KodjAj + A (djitj) + ^(A* - A 2 ) . (84) 

Time derivatives of tt q and A give rise to the following relations 

7T + = m*Ao , N 1 = . (85) 

Their time derivatives lead to Eqs. ([82]). The above four equations, however, are no 
longer enough to remove two extra degrees of freedom. Indeed, 

□ A + m*A° = -dun = 7T , (86) 

and since m 7 ^ 0, even if ttq were a plane wave solution of = 0, non-singular 
solutions A can be obtained. Therefore d^i is no longer constrained to zero, as it 
was in the massless case. Hence, such a theory propagates three degrees of freedom. 
Furthermore, using the relations (]K5j) one can rewrite the Hamiltonian as follows 



* ~ \«> + 5<**V + \< V 4 - " W ~ 4 + (87) 



Defining a new field 



we find the Hamiltonian 



D, = A, - ^, (88) 



where the pairs of canonical coordinates and their conjugate momenta are Bi and 
7Tj and A and The A field makes a negative contribution to the energy density. 
However, this field is decoupled from all the sources and, thus, cannot be produced 
to grow the negative energy. 
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Appendix B: Constrained non-Abelian gauge fields 



Here we generalize the constrained approach to the massless and massive non- 
Abelian fields. To understand the essence of these models we start with the discus- 
sion of the massless case. The Lagrangian reads as follows: 

I = F^,F atw + A^J^ + \ a {d^) a . (90) 
Equations of motion for the gauge fields and Lagrange multiplier are respectively: 

d v f^ + j; - d,x a = o , (91) 

d^A a » = . (92) 

An equation for A follows by taking a covariant derivative of ()91j1 

LTd^X" = D^Jl . (93) 

In a theory with a covariantly conserved source there could exist new non-trivial 
classical solutions for A such that <9 M A is also covariantly conserved. One particular 
solution is a path-ordered Wilson line 

d^X(x) = Vexp {-ig J A^z)dz^J M A(y) . (94) 

Thus, if there is a non-zero external colored current solution d^X at some point, then 
its value at any other point can be calculated according to 

The next issue to be address is that of quantum consistency of this approach. 
Again an important fact is that the Lagrangian (j50j) can be completed to a BRST 
invariant form: 

I = —-F^ v F afiV + A^J afl + A a (<9^) - ic a d»Dfc b . (95) 

Since at the classical solutions that we consider the FP ghost fields vanish, ()95|) 
recovers all the classical results of (|90|). The explicit BRST transformations leaving 
(J95|) invariant are: 

5Al = KDfc\ 5c a = - l -gCf abd c b c\ 5c a = \ a ( , 5X a = . (96) 

The expression for the BRST and Ghost currents, as well as the subsidiary conditions 
that guarantee the unitarity and completeness of the physical Hilbert space of states 
are the standard ones |4 . 

As in the case of a photon, the massless free propagator reads: 



Aff, = -5 ab Vtw d ^ dv,U . (07) 
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Finally, due to the presence of the BRST symmetry the Lagrange multiplier does 
not acquire kinetic term via the loop corrections. This is because the two-point 
correlation functions of the field is guarantied to be transverse due to the presence 
of the FP ghost and the transverse structure of the tree-level propagator (}9~Tj) 

cW (A* A h v )d v X h ~ & 1 X a 8 ah {r] lxv U - d li d v )d v X h = , (98) 

to all orders of perturbation theory 



B.l. Stiickelberg formalism 

Instead of using BRST symmetry to arrive from (|5U|l to ()95|). here we follow a 
conventional route. First we restore gauge symmetry of ()90|) and then fix that 
restored gauge invariance and introduce the appropriate FP ghosts. Let us start by 
defining new variables: 

igA li = U+DJJ, where U = e lt ^\ = «9 M + igB^ , (99) 

and rewrite the Lagrangian ()9(J|) in the following form 

, = U F ^ + ^HIES + E^lj,) . (10 o) 

2 V. 4 W W 



The above Lagrangian is gauge invariant under the local transformations of U — > 
e ia a t a jj anc | -q^ _^ e %a a t a g^ e -ia a t a _\_i^Q^a a t a j e -ia a t a ^ w j iere fa denote the generators 

of a local gauge group, and a a 's are gauge transformation functions 7 . We can proceed 
further and fix this gauge freedom by introducing into the Lagrangian the following 
term 

lTr(V^-£A) 2 . (101) 

This term removes quadratic mixing between and A, and after integrating out A 
and rewriting it back in terms of the A^ field, the result reads as follows: 



1 = l Tl [~\ F ^ Ftlv + y^ a ^ 2 + ^(V^X^) - icd^c + a^ 



(102) 



where 11^ = U A IM U + + ^[d^U]U + — A^. All the terms of the Lagrangian ()102j) . except 
the third one, is what one gets in the conventional approach. The third, gauge 
dependent term, has a structure that guarantees that the free propagator for an 
arbitrary value of £ coincides with the Landau gauge propagator of the conventional 
approach. This is due to the extra fields 7r a 's which by themselves have no kinetic 
term but acquire one through the gauge-parameter-dependent kinetic mixing term 
with A^ in (fTU2l . 

7 Notice that the current J M and the Lagrange multiplier A are not supposed to transform under 
these gauge transformations. This can be achieved by rewriting the fundamental fields out of which 
J M is constructed (as well as rewriting A) in terms of new fields rescaled by U's. Under the gauge 
transformations the new field transform in a conventional way but their variance is compensated 
by transformations of U's, so that and A stay invariant. 
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B.2. Comments on massive constrained non-Abelian fields 

Let us start with a component form of the constrained massive Lagrangian 

/ = ~F£ u F a ^ + A*J afl + \ a (d^) a + hvi 2 AlAl - ic a d»Dfc b . (103) 

It is interesting that the action of this theory is BRST invariant under the following 
transformations (the Lagrangian transforms as a total derivative) 

5A*=itDf£, 5c a = - l -gCf ahd c b c d , 5c a = \ a ( , 5\ a = iM 2 (c a . (104) 

Notice that the Lagrange multiplier does transforms w.r.t. the BRST and compen- 
sates for the non-invariance of the mass term. The above BRST transformations 
have the following peculiar properties: 

5 2 A£ = 5V = 0; 5 2 c a + + 5 2 X a ; 5 4 c a = 5 3 X a = 0. (105) 

The respective BRST and Ghost currents (the latter is due to the ghost-rescaling 
symmetry) c — > e a c, c — > e~ a c) can be derived 

jBRST = _ Fia/DuC + XE>imC _ llQ^fabdJ^ jGhost = i^ D ^a _ (^) C »)(106) 

The physical states could be defined in analogy with the standard approach 

Q B RS T |Phys >= 0, Q Ghost |Phys >= , (107) 

where the BRST and Ghost charges are defined as Qbrst = / d 3 xJQ RST (t, x) and 
Qchost = / d 3 xJQ host (t, x). This suggest that the A state, which upon the diagonal- 
ization of the Lagrangian acquires a ghost-like kinetic term, should not be a part of 
the physical Hilbert space of in and out states of the theory. The fictitious A par- 
ticle should be allowed to propagate as an intermediate state in Feynman diagrams 
softening the UV behavior of the theory, however, it cannot be emitted as a final 
in or out state of the theory. In this respect it should be similar to a FP ghost. 
However, because the peculiar properties of the above BRST transformations ()105|) . 
it still remains to be seen that for a rigorous construction of a positive semi-definite 
norm Hilbert space of states with unitary S-matrix elements the conditions (|107|) 
are enough. Besides the quantum effects, one should make sure that rapid classical 
instabilities are also removed. These may require further modification of the model. 
Detailed studies of this issue will be presented elsewhere. 

Having the formalism of the previous subsection developed it is easy now to 
discuss the spectrum of massive theory. First we restore the gauge invariance of 
the massive Lagrangian (|103|) by using the variables and then gauge fix it by 
(jl01|) . The resulting Lagrangian is 
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The mixing term between the gauge field and Goldstones is canceled and we integrate 
out the A field. The resulting Lagrangian reads: 

I = ^Tr ^F^(A) + l -M 2 Al + ^{d,A»Y + ^A^U,) + A,J^ ,(109) 

where, as before, we defined il M = UAJJ + + ^[d^U]U + - A^. 

A remarkable feature of this model is that the propagator takes the form: 

A ab = _ 5 ab ( _ (1+0 Q Q \ \ 5 ab frV^ (nQ) 

which has a smooth UV behavior and non-singular massless limit. 
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